Excited bands in odd-mass rare-earth nuclei 
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Abstract Normal parity bands are studied in 157 Gd, 163 Dy and 169 Tm using the pseudo SU(3) 
shell model. Energies and B(E2) transition strengths of states belonging to six low-lying rotational 
bands with the same parity in each nuclei are presented. The pseudo SU(3) basis includes states 
with pseudo-spin and 1 



and ^ and §, for even and odd number of nucleons, respectively. States 



with pseudo-spin 1 and | must be included for a proper description of some excited bands. 
PACS numbers: 21.60.Cs, 21.60.Fw, 23.20.Lv, 27.70.+q 



I. INTRODUCTION 

The nuclear shell model []IJ is the cornerstone in the 
microscopic study of nuclear properties. In the last two 
decades it has been possible to perform calculations for 
valence nucleons in the full sd— j^j and fp— || shells, 
providing a detailed description of energy levels, and elec- 
tromagnetic and weak transitions. Nevertheless, a shell 
model description of heavy nuclei requires further as- 
sumptions that include a systematic and proper trun- 
cation of the model space ||]. 

In light deformed nuclei the dominance of quadrupole- 
quadrupole interaction led to the introduction of the 
SU(3) shell model Q|. However, the strong spin-orbit 
interaction renders the SU(3) truncation scheme useless 
in heavier nuclei, while at the same time pseudo-spin 
emerges as a good symmetry || , and with it the pseudo- 
SU(3) model §. 

Pseudo-spin symmetry refers to the experimental fact 



that single-particle orbitals with j = I 



and j = 



(I — 2) + i in the shell r\ lie very close in energy, and can 
therefore be labeled as pseudo-spin doublets with quan- 
tum numbers j = j, rj = r\ — 1 , and 1 = 1 — 1. The origin 
of this symmetry has been traced back to the relativistic 
mean field equations ||. 

During the last decade the pseudo-SU(3) shell model 
has significantly evolved, becoming a powerful micro- 
scopic theory for the description of the low energy ro- 
tational bands in even-even heavy deformed nuclei, and 
of normal parity bands in heavy deformed odd-A nuclei. 
The first applications considered pseudo-SU(3) a dynam- 
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ical symmetry, using only one irreducible representation 
(irrep) of SU(3) to describe the yrast band up to the 
backbending regime [Q. After many improvements in 
the rotor Hamiltonian, a technical breakthrough came 
with the development of a computer code able to calcu- 
late reduced matrix elements of physical operators be- 
tween different SU(3) irreps M. From then on it was 
possible to use a realistic Hamiltonian which, besides the 
quadrupole-quadrupole interaction and rotor terms, in- 
cludes the breaking symmetry single particle energies and 
pairing interactions. 

In this way a fully microscopic description of low- 
energy bands in even-even and odd-A heavy deformed 
nuclei emerged. The model have been using as building 
blocks the pseudo-SU(3) proton and neutron states hav- 
ing pseudo-spin and i, which describe the even- and 
odd-A systems, respectively. The many-particle states 
are built as linear combinations of pseudo-SU(3) cou- 
pled states with well-defined particle number and to- 
tal angular momentum. Many rotational bands and 
B(E2) electromagnetic transition strengths have been de- 
scribed in the even-even rare earth isotopes 15e - 15& ^ 160 Qr\, 
i60,i62,i64 Dy and i64,i66,i68 Er g O, p) and in the odd 
mQOC 159,161 T >, i59,i63 Dy; ^E^^Tm, and 165 > 167 f, 



mass 
nuclei fll 



Tb 



Ei- 



Widening the landscape of applications of the pseudo 
SU(3) model, in the present work proton and neutron 
states having pseudo spin 1 and | , for even and odd num- 
ber of particles, respectively, are included in the Hilbert 
space. The final many-particle states in odd-mass nuclei 
have total spin |, | or S. In this enriched space, six or 
seven normal parity low-lying rotational bands in 157 Gd, 
163 Dy and 169 Tm are successfully described. Many of 
them have important pseudo-spin 1 and | components, 
evidencing their relevance in the present study. Intra 
and inter-band B(E2) transition strengths are also dis- 
cussed. This article complements the research of the scis- 
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sors mode, Ml transition strengths in the range between 
2 and 4 MeV, performed in the same nuclei JL7[ . 

In Section II a brief description of the pseudo SU(3) 
classification scheme is presented. The schematic Hamil- 
tonian and its parametrization is discussed in Section III. 
The results for low-lying rotational bands in 157 Gd, 163 Dy 
and 169 Tm are presented in Section IV, where their wave 
functions are also discussed. The B(E2) intra- and inter- 
band transition strengths are presented in Section V. The 
summary and conclusions are given in Section VI 



II. THE PSEUDO SU(3) BASIS 

The starting point in any application of the pseudo 
SU(3) model is the building of the many-body basis. The 
proton and neutron valence Nilsson single particle levels 
arc filled from below for a fixed deformation, allowing the 
determination of the occupancies in normal and unique 
parity orbitals |l5| . As it has been the case for all pseudo 
SU(3) studies up to now, nucleons in abnormal parity or- 
bital are frozen, being the dynamics described using only 
nucleons in normal parity states. This choice is further 
reflected through the use of effective charges to describe 
quadrupole electromagnetic transitions which are larger 
than those usually employed in shell- model calculations. 
While it has been shown that this is a reasonable and 
useful approach, it is nonetheless a strong assumption 
and the most important limitation of the present model. 
A more sophisticated treatment of the problem, with nu- 
cleons in intruder orbitals described in the same footing 
using SU(3) irreps is under development Jl8], |l9| . 

In Table | are presented the occupation numbers as- 
signed to each nuclei. 



the one with the lowest energy, is built as a linear com- 
bination 



Nuclei 


13 


rg 


n? 


r4 


nt 


157 Gd 


0.250 


8 


7 


6 


4 


163 Dy 


0.258 


10 


9 


6 


6 


169 Tm 


0.267 


11 


12 


8 


6 



TABLE I: Deformation and occupation numbers. 

The many-particle states of n a active nucleons in a 
given normal parity shell rj a , a — v or tt, can be classified 
by the following chains of groups: 

JV ~ - 

{!"<*} {fa} {f a }l a (A Q ,/i Q ) S 

U(n%) D U(Sl%/2) x U(2) D SU(3) x SU(2) D 
f T N 
50(3) x SU{2) D SUj(2), (1) 

where above each group the quantum numbers that char- 
acterize its irreps are given, and j a and n a are multiplic- 
ity labels of the indicated reductions. 

Any state |JjM), where J is the total angular mo- 
mentum, M its projection and i an integer index which 
enumerates the states with the same J, M starting from 



\JiM) =Y,C?\PJM) 
a 

of the strong coupled proton-neutron states 

\(3JM) 

\{U}(K^)S V , {f v }{K^ v )S v ; p(X(i)kL, S JM) 



(2) 



(LM L ,SM S \JM) (S-kM s „,S v M Sv \SMs) 



M L M S 

E 

k, v L v M v 



((A 77 £^77) ^77 -^77 



M v , S v Ms v ) 



(3) 



In the above expression (— ; — |— } and (— , — |— ) are the 
SU(3) and SU(2) Clebsch Gordan coefficients, respec- 
tively. 

In all previous applications of the pseudo SU(3) model, 
only those states with the highest spatial symmetry S ViV 
= and i have been included. In the present study states 
with iStt,!/ = 1 and | are also taken into account, allow- 
ing for coupled proton-neutron states with total pseudo- 
spin S = k, § or |. The pseudo-spin symmetry is close 
enough to an exact symmetry in atomic nuclei, support- 
ing a strong truncation of the Hilbert space. However, 
pseudo spin-orbit partners are not exactly degenerate, 
introducing a small degree of pseudo-spin mixing in the 
nuclear wave function. 

The quadrupole-quadrupole interaction can be ex- 
pressed in terms of the second order SU(3) Casimir op- 
erator C2, 



Q-Q = AC 2 - 3L{L + 1) 



(4) 



The eigenvalue of C2 for a given of SU(3) irrep (A,/i) is 
given by 

(C 2 ) = A 2 + fi 2 + Xfi + 3A + 3/i. (5) 

The larger the expectation value of C2, the greater the 
binding of that SU(3) irrep by a pure attractive Q ■ Q 
interaction. The pseudo SU(3) basis is built selecting 
those proton and neutron irreps with the largest (C2) and 
pseudo-spin and 1 , or | and 1 , for even and odd number 
of particles, respectively. The proton and neutron irreps 
are coupled to a total pseudo SU(3) (A, /x) irrep and to 
a total pseudo-spin S = |, | and |. The basis states 
employed in the present work for each of the three nuclei 
studied are presented in Tables |l[ III and IV. 



III. THE PSEUDO SU(3) HAMILTONIAN 

The Hamiltonian contains spherical Nilsson single- 
particle terms for protons {H sp 7r ) and neutrons {H sp u ), 
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total (A, fi) 




(A-7T, f-l^S-jr (A^, f-lv)Sv 


total (A, /i) 


(10,4)0 


(15,5)| 


(25,9)1 


(27,5)1 


(26,7)1 


(7,7)1 


(24,0)0 


(31,7)1 








(10,4)0 


(18,2)1 


(27,5)1 


(26,7)1 


(29,4)1 


(H,2)i 


(20,5)0,1 


(31,7)1,| 


(33,3)1 


3 
2 


(32,5) 2 ,§ 






(28,6)1 


(30,2)1 




(7,7) \ 


(16,10)0 


(23,17)1 








(10,4)0 


(16,3)|,| 


(27,5)1 


2 (26,7)^,2 




(H,2)| 


(24,0)0 


(33,3)1 


(34,1)1 




(35,2)1 


(10,4)0 


(17,1)|,| 


(27,5)1 


3 
2 




(11,2)| 


(21,3)0,1 


(33,3) 2 , 2 


(32,5)1 


3 
2 


(34,l)i,f 


(11,2)1 


(18,2)| 


(27,5)1 


2 (29,4)^,2 




(H,2)| 


(22,1)1 


(33,3) 2 1 2 


(34,1)1 


3 
2 








(31,0)1 


2 (30,2)2,2 




(9,3) 2,5 


(24,0)0 


(33,3) 252 








(11,2)1 


(15,5)| 


(27,5)1 


2 (26,7)^,2 




(8,5)i,| 


(24,0)0 


(32,5) 2 1 f 








(11,2)1 


(16,3)1, § 


(27,5)1 


3 5 
2 ' 2 




(io,i)| 


(24,0)0 


(34,1)1 








(12,0)0 


(15,5)| 


(27,5)1 






(7,7)| 


(20,5)0,1 


(27,12)1,- 


i 
> 






(12,0)0 


(18,2)| 


(30,2)1 






(H,2)| 


(16,10)0 


(27,12)1 








(8,5)0,1 


(18,2)| 


(26,7)1 


I (27,5)2,| 




(4,10)1 


(16,10)0 


(20,20)1 








(9,3)0,1 


(18,2)1 


(27,5)1 


3 
2 

















TABLE II: The pseudo SU(3) irreps used in the description 
of 157 Gd. 





total (A, /i) 


(10,4)0 


(16,7)1 


(26,11)1 






(7,7)0,1 


(19,4)1 


(26,ll)i,| 






(10,4)0 


(19,4)1 


(29,8)1 


(31,4)1 


(30,6)1 






(32,2)1 






(10,4)0 


(20,2)1 


(31,4)1 


(30,6)1 


(32,2)1 


(10,4)0 


(21,0)| 


(31,4)| 






(11,2)1 


(19,4)1 


(31,4)i,| 


(30,6)i,| 


(32,2)i,| 


(11,2)1 


(20,2)1 


(31,4)i,| 


(32,2)1,| 




(12,0)0 


(19,4)| 


(31,4)1 






(11,2)1 


(21,0)| 


(32,2)1,1, 


5 
2 




(12,0)0 


(20,2)1 


(32,2)1 







TABLE III: The pseudo SU(3) irreps used in the description 
of 163 Dy. 



the quadrupole-quadrupole (Q ■ Q) and pairing interac- 
tions (Hp a i r ^ and H va ir,v)i as well as three 'rotor- like' 
terms which are diagonal in the SU(3) basis. 



H — {Hsp.u — Ga H pa i r}a } — - \ Q ' Q (6) 

a — 1Z ,V 

+ a Kj + b J 2 + A sym C 2 . 

This Hamiltonian can be separated into two parts: the 
first row includes Nilsson single-particle energies and the 
pairing and quadrupole-quadrupole interactions (Q is the 
quadrupole operator in the pseudo SU(3) space, see be- 
low). They are the basic components of any realistic 
Hamiltonian and have been widely studied in the nuclear 
physics literature, allowing their respective strengths to 
be fixed by systematics [^0| The SU(3) mixing is due 
to the single-particle and pairing terms. The second row 
of the Hamiltonian (IfJ) contains the so-called 'rotor-like' 



TABLE IV: The pseudo SU(3) irreps used in the description 
of 169 Tm. 



Nuclei 


a 


b 




157 Gd 


0.046 


0.0020 


0.0008 


163 Dy 


-0.040 


0.0040 


0.0016 


169 Tm 


0.019 


0.0007 






TABLE V: Parameters of the Hamiltonian 



terms, used to fine tune the moment of inertia and the 
position of the different K bands. They have been stud- 
ied in detail in previous papers where the pseudo SU(3) 
symmetry was used as a dynamical symmetry Q. The 
strength of these three terms were the only ones adjusted 
nuclei by nuclei. A detailed analysis of each term of this 
Hamiltonian and its parametrization has been presented 
elsewhere In Table |v| are shown the current values. 
These are the same used before for these three nuclei [^2| , 
when the proton and neutron subspaces were reduced to 
S = and \ (in other words, the same set of parameters 
are working in both approximations of the theory) . 
The electric quadrupole operator is expressed as |?J 



(7) 



with effective charges = 2.3, e„ = 1.3. These values 
are very similar to those used in the pseudo SU(3) de- 
scription of even-even nuclei 0, ^3| . They are larger than 
those used in standard calculations of B(E2) strengths 
p0| due to the passive role assigned to the nucleons 
in unique parity orbitals, whose contribution to the 
quadrupole moments is parametrized in this way. 

The inclusion of configurations with pseudo-spin 1 and 
I in the Hilbert space allows for a description of several 
highly excited rotational bands in odd-mass nuclei. Their 
effect on the Ml transition strengths have been discussed 
in Ref. |rij| . The most relevant result reported there is 
that, when the configuration space is restricted to states 
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with pseudo-spin and \, it is not possible to find p^] 
any Ml excitations with sizable strength from the ground 
state toward the states between 2 to 4 MeV energy region 
for any of the three nuclei. The Ml strength appears only 
when the Hilbert space is enlarged to include states with 
S„ or S v = 1 or 3/2. 



IV. ROTATIONAL BANDS AND SPIN 
CONTENT 

Fig. |l| shows the yrast and excited normal parity bands 
in 157 Gd. Experimental p3 data are plotted on the left 
hand size of each column, while those obtained using the 
Hilbert space and the Hamiltonian parameters discussed 
in the previous sections are shown in the right hand side. 
The agreement between both is excellent, but for the two 
higher energy bands the lack of more experimental data 
in bands E and F (see fig. [I]) prevents a more rigorous 
comparison. 
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FIG. 1: Negative parity bands in 157 Gd. The integer numbers 
indicate twice the total angular momentum of each level. 



The whole energy spectra is built up by the interplay 
between the single-particle and quadrupole-quadrupole 
terms in the Hamiltonian. The use of realistic single- 
particle energies plays a key role in the appropriate or- 
dering of the different band-heads. 

Figs. |^ and |^ show the normal parity bands in 163 Dy 
and 169 Tm, respectively. In 163 Dy the agreement be- 
tween theory and experiment is very good for the seven 
rotational bands (A-G). These seven bands represents 
nearly all the bands measured. For 169 Tm the agree- 
ment between the theoretical band structure and its ex- 
perimental counterpart is still good, although the differ- 



ences in bandhead energies are larger than in the previous 
cases. 
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FIG. 2: Negative parity bands in 163 Dy. 
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FIG. 3: Positive parity bands in 169 Tm. 

The bands plotted in Figs. (|lH|) have a very regu- 
lar structure as one moves up the bands. The mixing- 
remains nearly the same for the states with different an- 
gular moments belonging to the same band. As a con- 
sequence, the pseudo-spin content is practically constant 
along all the members of each band. Table |Vl] shows the 
calculated pseudo-spin content of each rotational band 
in 157 Gd, 163 Dy and 169 Tm, respectively. As shown in 
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Nuclei 


Band 


S = 1/2 (%) 


S = 3/2 (%) 


157 Gd 


A 


89 


11 




B 


92 


8 




C 


92 


8 




D 


3 


97 




* E 


77 


23 




F 


16 


84 


163 Dy 


A 


100 







B 


100 







C 


100 







D 


100 







E 





100 




* p 


100 







G 





100 


169 Tm 


A 


93 


7 




B 


100 







C 


9 


91 




D 


100 







E 


35 


65 




F 


37 


63 



TABLE VI: Pseudo-spin content for each band in Gd, 
163 Dy and 169 Tm. In first column it is indicated the nuclei, 
the second shows the band, and third and four columns show 
the 1/2 and 3/2 pseudo spin-content as percentage for each 
band, respectively. 
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FIG. 4: B(E2) intra-band transitions in 157 Gd [e 2 b 2 ]. 



Eq. 



each eigenstate is described as a linear combina- 



tion of pseudo SU(3) states \/3JM), with pseudo-spin |, 
| and |, depending on the proton and neutron pseudo- 
spin couplings. 

In most of the bands the mixing is very small, about 
half of them have no mixing at all. The main contribution 
of the present work is to add those states with pseudo- 
spin | (and those marked with a *, see below) to the 
pseudo SU(3) description of these nuclei. 

The total pseudo-spin content of the nuclear wave func- 
tion is built with the coupling of the SV and S„ compo- 
nents. In all the cases, the states with pseudo-spin S = I 
are built as 1 ® h, and the ® | is not present in the 
bands analysed. By the other side, although most of the 
S = h come from the coupling 

with a * in Table VI, are built as SV = 1 



|, the bands marked 
i for 



S u 



the band E in 157 Gd and = | (g> S u = 1 for the band 



mi 



Dy). 



V. B(E2) TRANSITION STRENGTHS 

In Figures ||, || and ^ are shown the B(E2) intra-band 
transition strengths for the reported normal parity bands 



in 57 Gd, bJ Dy and w Tm, respectively. All experimen- 
tal data available |24j have been included with error bars 
for comparison. For a given band, labeled A to G, the x 
axis represents twice the angular momentum of the ini- 
tial state for each transition, while the y axis shows the 
B(E2) strength. The graphs on the left hand side show 
B(E2;J — » J+l) transitions, and those in the right hand 
side present the B(E2;J — » J+2), both in units of e 2 b 2 . 
The "+" symbols with segmented lines refers to the cal- 
culated values, and the diamonds with error bars to the 
experimental ones. 



In Tables VII , VIII and [X are shown many inter-band 
B(E2) transition strengths in 157 Gd, 163 Dy and 169 Tm, 
respectively. All known experimental data available 
were included between parenthesis. Theoretical values 
are in agreement with those reported previously for the 
lowest energy bands in 163 Dy [ fhS] . 

The agreement between the calculated values with 
those measured is remarkable. Some transitions, like 
3/2" -> 7/2-, 3/2 D -> 7/2- in 157 Gd, and 3/2 B -> 
7/2", 3/2^ -» 7/2", 9/2 c -> ll/2 B in 163 Dy, are larger 
than 0.05 - 0.10 e 2 b 2 . It reflects the mixing between the 
different bands. 

The B(M1;J£ S -> J J) transitions in 157 Gd, 163 Dy and 
169 Tm have been presented in Ref. p~~ 
the ground states | , 



and i + , respectively, in these 



Jg s refers to 
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TABLE VII: Calculated B(E2) inter-band transitions in 
157 Gd [eV x 10- 2 ]. 
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TABLE VIII: B(E2) inter-band transitions in 163 Dy. The 
experimental data are given between parenthesis [e 2 b 2 x 10 -2 ]. 
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FIG. 5: B(E2) intra-band transitions in lea Dy [e V 



nuclei. In that work it was shown that most of the states 
with energies between 2 and 4 MeV have very impor- 
tant contributions from states with proton and neutron 



pseudo-spin 1 and 



VI. SUMMARY AND CONCLUSIONS 

The pseudo SU(3) shell model for odd-mass nuclei has 
been shown to offer a quantitative microscopic descrip- 
tion of several normal parity rotational bands in 157 Gd, 
163 Dy and 169 Tm. The present article complements the 
study of the scissors mode in these nuclei, their fragmen- 
tation and their summed B(Ml;j) strengths. 

In order to successfully describe the excited bands, 
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FIG. 6: B(E2) intra-band transitions in 169 Tm [e 2 b 2 



intra- and inter-band B(E2) and B(M1) transition 
strengths, it was necessary to enlarge the Hilbert space, 
including those pseudo SU(3) irreps with the largest C2 
values having pseudo-spin 1 and | and to use realistic val- 
ues for the single-particle energies. This expansion of the 
model space allowed for the description of new low lying 
normal parity bands, in agreement with the experimental 
data. The new bands have predominantly S = I, but it 
was shown that the most important contribution comes 
from the proton or neutron subspaces with S = 1. It 
implies that the pseudo-spin mixing in the wave function 
takes place mostly in the sub-space with even number of 
particles. The interplay between the single-particle and 
the quadrupole-quadrupole terms in the Hamiltonian de- 
fines this mixing. The delicated balance between these 



two interactions, defines the gross features of the calcu- 
lated excited bands. 

Intra- and inter-band B(E2) transition strengths are 
predicted in good agreement with their measured val- 
ues and with those reported previously with the pseudo 
SU(3) model. Detailed description of magnetic proper- 
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TABLE IX: B(E2) inter-band transitions in 169 Tm. Experi- 
mental data available are between parenthesis [e 2 b 2 x 10 -2 ] 



ties of these rotational bands |25| and the excited bands 
in other rare-earth nuclei are desirable. 
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